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Abstract, a clustering and a condensation in lighter nuclei is presently strongly and 
increasingly discussed in the literature both from the experimental side as from the theoretical 
one. A discussion of the present status of the theory as well as outlooks for future developements 
will be presented. 



1. Introduction 

Since about ten years when the idea of a possible existence of a condensate type of states in na 
nuclei was advanced and formulated for the first time many exciting new results, theoretical 
and experimental ones, have been produced. In this contribution, we would like to assess where 
we stand and what possible future extensions may be. 

Let us start with the reminder that nuclear clustering and in particular a clustering would 
not exist, if we did not have in nuclear physics four different types of fermions (proton/neutron 
spin up/down), all attracting one another. We should be aware of the fact that this is a rather 
singular situation in fermionic many body systems. However, the possibility of future trapping 
of four different kinds of cold fermionic atoms, may open a new field of cluster physics with 
similar features as in nuclei. In a mean field description of an isolated a particle (what with, 
e.g., Skyrme forces gives reasonable results, if the c.o.m. motion is treated correctly) the four 
fermions can occupy the lowest OS-level. Would there be only neutrons, only two of them can 
be in the OS-level, the other two neutrons would have to be in the energetically very penalising 
OP-state. That is why a particles exist, tetra- neutrons not. The ensuing fact is that a particles 
are very strongly bound (E/A ~ 7 MeV), almost as strong as the most bound nucleus which is 
56 Fe (E/A ~ 8 MeV). In addition the first excited state of the a particle (~ 20 MeV) is by factors 
higher than that of any other nucleus. The a particle can, therefore, be considered as an almost 
inert ideal bosonic particle. As we will see in the discussion below, inspite of its strong binding, 
a particle condensation can only exist in the so-called BEC (Bose-Einstein Condensation) phase 
what implies low density. There is no analogue to the BCS phase of pairing where the Cooper 
pairs can have very large extensions, strongly overlapping with one another, still being fully 
antisymmetrised. This is the reason why a condensation only can be present at low densities 
where the a particles do not overlap strongly (this holds, if the system consists of protons and 
neutrons and q's. If other clusters as t, 3 He,d are around, the situation may change, see below). 



These considerations apply for nuclear matter as well as for finite nuclei. The Hoyle state in C 
which can to a good approximation be described as a product of three a particles occupying all 
the lowest OS state of their bosonic mean field has a density which is by a factor 3-4 lower than 
the one of the ground state of 12 C. In the ground state there exist a- type of correlations but 
there is no condensation phenomenon. Let us start our considerations with infinite matter. 

2. a particle condensation in infinite matter. 

The in medium four-body equation can be written in the following form 

(E a , K -£ 1 -E 2 -£ 3 - £4)^1234 = (1 - /l - /2)wi21'2'*i/ 2 /34 + (1 - /l - /3)ui31'3'*i/ 23 '4 + - W 

In total, there are six terms coming from permutations. The £j are kinetic energies plus mean 
field corrections; ^1234 are the matrixelements of the two body interaction, and fa is a Fermi- 
Dirac distribution of the uncorrelated nucleons accounting for phase space blocking. Repeated 
indices are summed over and index numbers comprise momenta and spins. The above equation 
considers, therefore, one quartet in a gas of uncorrelated nucleons at temperature T. The 
analogous two body equation can be used to determine the critical temperature T c for the onset 
of superfluidity or supraconductivity where T c has to be determined so that the eigenvalue comes 
at two times the chemical potential /i. This is the famous Thouless criterion of BCS theory. In 
analogy with pairing, one has to find the critical temperature T C Q so that the eigenvalue of the 
four body equation (1) comes at 4/i. The in medium four body equation is very difficult to solve. 
Nontheless, the solution has been found employing the Faddeev-Yakubovsky equations and using 
the Malfliet-Tjohn bare nucleon-nucleon interaction which yields realistic nucleon-nucleon phase 
shifts and properties of an isolated a particle [2]. To simplify the problem, we made in addition 
a very easy to handle variational ansatz of the four body wave function in (1). It consists of a 
mean field ansatz for the a particle projected on good total momentum. In momentum space 
this is 

1-1234 oc S(K - ki - k 2 - k 3 - k4)^(k 1 )v?(k 2 )^(k 3 )99(k4) (2) 

Inserting this ansatz into (1), one obtains a nonlinear HF-type of equation for the S-wave 
function </?(k). Of course, for quartet condensation, we choose K = 0. With the mean field 
ansatz (2), one cannot use a bare force. We adjusted an effective separable force with two 
parameters which are chosen to reproduce binding energy and radius of the free a particle. 
The full Faddeev-Yakubovsky solution of (1) is shown for symmetric and asymmetric matter 
in Fig. 1 (crosses). We see that the ansatz (2) which very much eases the otherwise difficult 
solution of (1) works very well (continuous line). Also shown is the critical temperature for 
deuteron condensation. The striking feature is that a particle condensation abruptly breaks 
down already at very low density which approximately coincides with the point where the a's 
start to overlap appreciately (this fact was already found in [3] using a somewhat different 
variational ansatz for ^1234). On the other hand deuteron condensation goes on up to very 
high densities and the limit is only triggered by the range of the effective force (which has been 
readjusted to reproduce the deuteron properties). This is so for symmetric nuclear matter. For 
strong asymmetry, deuteron condensation breaks down earlier than a condensation because the 
a particle according to its much stronger binding is less sensitive to asymmetry This is the 
afore mentioned phenomenon that a condensation only exists in the BEC phase, i.e., at low 
density, whereas deuteron continuously goes from negative to positive chemical potentials where 
for the latter the deuterons turn into large size Cooper pairs. More on this can be found in [2] 
[3] . We should mention that our calculation of T" is only reliable rather close to the break down 
point. For lower densities, the should join the one for condensation of ideal bosons (a's). 
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Figure 1. Critical temperatures for a 
particle and deuteron condensation in 
symmetric nuclear matter (IS (I function 
of \i (a) and density (b). 



To describe this feature, one should extend our theory to the so-called Nozieres Schmitt-Rink 
(NSR) theory [5] for pairing, see also [5], to a particle condensation. This, however, has not 
been worked out so far and remains a task for the future. 

At zero temperature, there are many a's which go into the condensate phase. For this, we 
have to set up an approach analogous to the nonlinear BCS theory. Equation (1) corresponds 
to the linearised version and only describes one a particle in an otherwise uncorrelated gas (at 
finite T) of fermions. In finite nuclei, there may exist such a situation even at zero temperature. 
This is the case of 212 Po which can to a certain extent be viewed as an a particle sitting on top 
of the doubly magic core of 208 Pb which can be well described by a HF-mean field approach, 
i.e. a Fermi gas in a container. We will come back to this later when we discuss finite nuclei. 

After having considered the linearised version of the equation for the quartet order parameter 
at the critical temperature, let us now try to write down, in analogy to the BCS case, the fully 
non-linear system of equations for the quartet order parameter. To clearly see the analogy to 
the BCS case, let us repeat the latter equations in a slightly unusual way. The pairing order 
parameter allowing for non-zero c.o.m. momentum of the pairs, reads 



(e/ci + £fc 2 ) K kik 2 + (1 - n kl - n fc2 )A klk2 = 2 / uAt klk2 (3) 

where A klka = £ vu^h^, «k i; k 2 , and K klk2 = (BCS|c kl c ka |BCS) = u kl f ka is the usual 
pairing tensor (with spin indices suppressed) and n k = v\ = 1 — u\ are the BCS occupation 
numbers. The e^'s are, as before, the kinetic energies, eventually including a mean field 
correction. The occupation numbers can be obtained from the Dyson equation 

Gt = G° k + G° k M%G% with M k = ^± (4) 

the BCS mass operator where A/% is the diagonal part of the gap for cases where the pairs 
are at rest. From the single particle Green's function, obviously we can calculate the occupation 
numbers closing, thus, the typical BCS selfconsistency cycle. Inspired by the BCS case, we then 
write for the quartet order parameter, see (1) 



(4/i - El - £ 2 - £3 - £4)K1234 = (1 - Til - n 2 )t>i21'2'K4'2'34 + (1 - Til - n3)ui31'3'Kl'23'4 + ••• (5) 



with K4234 = (C1C2C3C4}. Again, we have to close the equation with the Dyson equation for 
the occupation numbers. However, the massoperator now contains the quartet order parameter 



= y A 1234 [n^o^ + nOnK]At 234 

^34 w + £2 + £3 + £4 

with n° = 1 — n° and being the uncorrelated Fermi function, i.e., the Fermi step and 
A1234 = Yl ^i23'4'^3 / 4'34 where V1234 is an effective coupling vertex linking the single particle 
motion to the order parameter (for more details and derivation, see [7]). Before trying to solve 
this equation, let us discuss the differences between the pairing and the quartet case. The first 
thing which strikes is that the three 'holes' only have to have total momentum k2+k3+k 4 = — ki 
and, therefore, we have a remaining sum over momenta. In the pairing case with only one 'hole', 
there is no sum. Furthermore in the pairing case the hole propagator has no phase space factor 
because 'forward' and 'backward' going parts add up to one: n° + n° = 1. In the quartet case 
there are three hole propagators and the corresponding sum of phase space factors does NOT 
add up to one, i.e. n^n^n^ + n^n^n® 7^ 1!! This makes a dramatic difference with the pairing 
case. In order to understand this a little better, let us compare the level density of a single hole 
with the one of three holes: 

g lh (oj) = J>° + n° k ]6(u + £ fc ) = E ^ + £ fc) ( 7 ) 
k k 

g 3h (u)= K« + <«] ,5 ( w + e fei + ^ 2 + e fc3) (8) 

kik 2 k 3 

The 3h level density is shown in Fig. 2. We see that for positive \x there is a striking difference 




Figure 2. 3h level densi- 
ties for negative and positive 
chemical potentials, respec- 
tively. Note that on the hor- 
izontal axis, the origine cor- 
responds to fj, = 0. 



with the lh level density. At positive /x, g lh (u) = fi) is obviously finite (not shown), whereas 
g 3h (uj = 3/u) goes through zero. This is because phase space constraints and energy conservation 
cannot be fullfilled simultaneously at the Fermi surface in the latter case as is easily verified. 
That is, in the quartet case, exactly at the point where the correlations should be built up, 
namely at the Fermi-level, there is no level density! As a consequence, no quartet condensation 
is possible for positive /u. On the contrary for negative fj,, n° = and, thus, the phase space 
factor in the case of g 3h is also equal to one and then there is no qualitative difference with the 
lh case. This explains in a natural way why quartet condensation is not possible at positive 
chemical potential. It is by the way well known that any mp — nh level density, besides the 
single particle case, goes through zero at the Fermi level. For example, the origin of fermions at 
the Fermi energy having an infinite free mean path stems from the fact that the 2p-lh (2h-lp) 
level density is zero at the Fermi level what is the equivalent to the imaginary part of the optical 
potential being zero at that energy. Also the 2p-2h level density which plays an important 
role for the damping of zero sound modes is zero at zero frequency. As first pointed out by 
Landau, it approaches zero as uj 2 . In conclusion for positive [i only pairing survives whereas 



quartetting breaks down and only exists in the BEC phase with negative p. The situation may 
be different when other light clusters are present, i.e. in a mixed gas of, e.g., nucleons, tritons 
( 3 He), deuterons. Then a nucleon with momentum k may, eventually, directly pair with, e.g., 
a triton of momentum — k (or the other way round), rather similar to the standard pairing 
situation besides the fact that now two fermions with different masses have to pair up. Similar 
considerations hold for the pairing of two deuterons (pairing of 'bosons'). In compact star 
physics such situations may exist when the star is cooling down. The extension of our theory to 
this scenario is a task for the future. The full solution of the nonlinear set of equations (6) and 
(7) is again very much eased in taking for the order parameter the factorisation ansatz (2). The 
most interesting result is that the occupation numbers, e.g., for p around zero is far from being 
close to saturation. At k = it is only approximately n^. = o ~ 0.30. This scenario is analogous 
to pairing in the BEC regime. More results can be found in [7j. 



3. Finite Nuclei 

As we know from pairing, a direct observation of condensation phenomena only is possible 
in finite nuclei. Of course, in such small systems, there cannot exist a condensation in the 
macroscopic sense. Nevertheless, as we know very well, only a handfull of Cooper pairs suffices 
to show clear signatures of pairing in nuclei. For a particle condensation it is the same story. 
We also only can expect that there exists about a handfull of a particles, essentially in lighter 
na nuclei, in a gaseous phase at low density. It is, indeed, surprising that such states at low 
density with about p = po/3 — po/4 with po the density at saturation, indeed exist as excited 
quite long lived states in those nuclei. The most fameous example is the Hoyle state in 12 C at 
7.65 MeV, just about 300 keV above the 3a threshold. We will not dwell much on the successful 
theoretical description of this state (and others,e.g., in 16 0) with the THSR wave function [T], 
since this has been presented in the recent literature a great number of times [8]. Let us only 
make a couple of remarks. The THSR wave function is schematically written for a finite number 
of quartets as 

*™SR = A [$ a $ a _$ a ] (9) 

where the single a wave function <3? Q depends on four spatial coordinates, the spin-isospin 
part being suppressed. This wave function is fully antisymmetric due to the antisymmetriser A 
and is analogous to the number projected BCS wave function 

* BCS = ^ pair </> pair pair ] (10) 

where 4>p a [ T is the Cooper pair wave function depending on two spatial coordinates. The 
calculus with the a condensate wave function is very much facilitated with a variational ansatz 
where <& Q is split into a product of a c.o.m. gaussian with a large width parameter B times 
another, intrinsic, gaussian depending only on the relative coordinates of the a particle and 
having a width parameter b which corresponds to the size of an isolated a particle. The first 
important remark to be made is that this THSR wave function contains two important limits: 
if B = b then it corresponds to a pure harmonic oscillator Slater determinant. If B » b, then 
the a's are so distant from one another that the Pauli principle among the different a's can be 
neglected and, thus, the antisymmetriser be dropped. The THSR wave function is then a pure 
product state of a particles, i.e., a condensate state of ideal bosons. Reality is, of course, in 
between those limits and one important task is to find out whether reality is closer to a Slater 
determinant or to a Bose condensate. That this question must be carefully investigated, can 
also be deduced from the fact that a number projected BCS wave function always leads to a 
non trivial pairing solution. For example for 208 Pb, one obtains a non trivial BCS solution 
inspite of the fact that 208 Pb is certainly not superfluid. Only when the original particle number 



breaking BCS theory with |BCS) = e^ Zkk ' CkC - k \vac) has a non-trivial solution, we can speek 
of a superfluid nucleus. For 208 Pb there is no such solution. One way to analyse whether the 
THSR approach leads primarily to an a condensate or to a Slater determinant, is to investigate 
the bosonic occupation numbers. An ideal Fermi gas has occupation numbers which are one or 
zero. In an ideal Bose condensate the bosons will occupy the lowest single particle state with 
100 %. Of course, in real nuclei, neither the fermionic occupation numbers nor the bosonic 
occupations are the ones of an ideal gas. For nucleons, we know that occupation numbers are 
around 70-80% down from unity. For the a occupations in the condensate, we also have found 
a number around 70%, all the other occupancies being down by at least a factor of ten, whereas 
for the ground state the occupation numbers are almost equally distributed according to the 
SU3 shell model scheme [8] . This is a typical condensate situation, though not totally the 
one of an ideal condensate. Residual antisymmetry effects between the slightly overlapping a 
particles makes that with a certain probability the a's are scattered out of the condensate. In 
addition there are also other correlations at work. One of the most important one is certainly 
the formation of 8 Be clusters out of two a's. It is even not clear whether an a gas does not in 
reality consist out of a gas of 8 Be's. The latter are, of course also bosons and the old question 
arises whether in an attractive Bose gas the bosons condensate as singles or as molecules |9]. 
This is certainly a very interesting question which desserves further studies in the future. 

Many other extensions of a condensation are presently discussed theoretically and 
experimentally. An intersting aspect is whether on top of the a condensate states excited a 
gas states exist. A long debate has recently been closed about the nature and existence of the 
second 2 + state close to 10 MeV excitation energy in 12 C. Very nice experimental results by M. 
Freer, M. Itoh, and M. Gai have recently shown that this state is there and that it is part of 
a family of a gas states |10||llj. Many more results, for instance, in 16 are to be expected. 
One of the most exciting aspects is that one may be able to dismantle na nuclei with n a rather 
large number like n = 10 or more into n a particles. First results in this direction have been 
reported at this conference with 56 Ni, i.e. a nucleus with 14 a particles by H. Akimune. A dream 
would be that all a particles be just excited to the Ikeda threshold and then they disintegrate 
in a very slow motion as a kind of coherent state driven by the Coulomb force, i.e., it would 
be some kind of soft Coulomb explosion. This would then be rather close to what happens 
with a trapped Bose condensate of cold atoms upon switching off the trapping potential. Other 
exciting perspectives are that a particles could exist in a gaseous phase on top of an inert core. 
For example four a's on top of 16 in 32 S, or other variants, even in quite heavy nuclei. Because 
of space restrictions, we, unfortunately, cannot extend further on these exciting aspects of a gas 
type of states in nuclei. 

Before closing, we would like to discuss, however, the question of a possible preformation of a 
particles in heavy nuclei. As well known, this question is of great importance fo the description 
of a decay rates. Let us take the example already alluded to in the infinite matter section, 
namely 212 Po. Since the lead core is doubly magic, one can treat it in mean field, i.e. as a Slater 
determinant with a suitable Skyrme or Gogny type of force. That is we may view 212 Po as an 
a particle on top of a finite Fermi gas. If such a configuration exists, and the a decay rate of 
212 Po suggests that, it is clear from our experience with the Hoyle state that such a cluster state 
cannot be described within the shell model alone, see [13] [H]. In the infinite matter section, 
we have already considered a situation where a single a particle is imbedded in an uncorrelated 
Fermi gas. This was the case with respect to the critical temperature. For 212 Po we can consider 
our treatment of the single a case even at zero temperature. Of course we cannot use the ansatz 
(2) as is, since it corresponds to a c.o.m. wave function which is a plane wave e* KR . In position 
space (2) namely reads ^1234 oc e tKB "ipi n t.(\ r i — r j\) where ipi n t. is the intrinsic wave function 
depending only on the relative coordinates. In a finite nucleus, instead of the plane wave, one 
would have to use a wave packet considering the following wave function 



^1234 oc $(R)^ nt .(|ri - rj \) (11) 

where ipi n t. and <3? are to be determined variationally upon inserting (12) into eq (1) where all 
indices and ingredients correspond now to the mean field of 208 Pb. For instance the £j correspond 
to the HF energies and the fi are to be replaced by the HF occupation numbers at T = 0. In 
addition for the ifiint. , one could use, as for the study of the Hoyle state, a gaussian wave function 
with the width parameter b. Then the single unknown would be the spherical wave function 
&(R) to be determined variationally where R is the c.o.m. coordinate of the a with respect to 
the center of the Pb core. Our knowledge that a's only can exist at very low density incites 
us to believe that this $ wave function should be peaked rather far out in the surface of the 
Pb nucleus. If true, this would be a nice explanation of a preformed a particle in the nuclear 
surface. Adding more than one a to the Pb core may suggest that there could exist some sort 
of a condensate in the surface in a fluctuating state. However, as we know, adding more a's to 
the Pb core leads to deformed nuclei and the treatment of a's in the surface becomes a much 
more delicate subject. 

4. Conclusions 

As we have seen, the existence of a gas and a condensate states in nuclear systems where the a's 
play practically the role of elementary bosons, is fascinating. Nuclear physics is at the forefront 
of this kind of physics. In the future, experiments with cold atoms trapping four (or more) 
different kinds of fermions may also open wide perspectives in the field of cluster physics. For 
more reading on a cluster states, we invite the reader to consult our review article [8]. 
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